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Abstract 

The charmonium transitions from the XcJ states to rj c with emission of one or 
two pions are considered. It is shown that the only processes of such type arising 
in the leading order in the QCD multipole expansion are the decays Xd i] c 7r + ir~ 
and XcO — > Vc^ - The absolute rate of the latter decay, in spite of being isotopically 
suppressed, is by about an order of magnitude larger than that of the former, so that 
the branching fractions for the two processes are approximately equal, and are expected 
to be significantly below the recent upper bound reported by BESIII. The rates of both 
decays are also related to that of another similar transition h c —¥ J/tp-rr . 



Recently the BESIII experiment has established [T] interesting upper bounds on the rates 
of two-pion transitions from the Xcj charmonium states to the spin-singlet ground state r] c , 
Xcj r] c 7T + 'n-~ . The hadronic transitions between states of heavy quarkonium with emission 
of a pion pair as well as of isospin-violating transitions with emission of a single neutral 
pion attract a considerable and long-lasting interest as being related to understanding the 
QCD dynamics of heavy and light hadrons. The purpose of the present paper is to revisit 
the theoretical description of such processes within the multipole expansion in QCD (2J [3], 
specifically in application to the transitions between P-wave and .S-wave states of a heavy 
quarkonium, which may be helpful in assessing the possibilities for further experimental 
studies of such transitions. The main result to be argued here is that out of four allowed 
transitions 0, Xc2 Tj c n + ^~ ', Xc2 Vc^ Xd ~~ >* 7 7c 7r+7r_ ) an d Xco — > Vc^ , only the latter two 
are contributed by the leading El — Ml term of the QCD multipole expansion, so that both 
transitions from the J = 2 state Xc2 should be strongly suppressed. Moreover, the absolute 
rate of the isospin breaking process Xco — >■ Vc^ should be larger than that for the decay 
Xci — >■ r] c 7T + 7r~ by approximately an order of magnitude, which essentially compensates for 
similar difference of the total widths between Xco and x& an d makes the branching fractions 
for the two transitions approximately equal: 



Furthermore, the absolute rate of the transition Xco — > Vc^ is simply related to that for a 
similar decay of the spin-singlet P-wave charmonium h c : 



One more remaining transition of the same origin, the decay h c — > J/ipirn, is strongly 
suppressed kinematically, and is unlikely to be of an immediate significance for experimental 
studies. 

Neither of the discussed IP — > IS transitions in charmonium has been observed so 
far. (With a possible exception of h c — > J/ifm°, whose sighting reported by the E760 
experiment [5] has not been confirmed.) A theoretical expectation for the absolute values 
of the discussed decay rates is uncertain due to poor knowledge of the relevant quarkonium 
transition amplitude for the IP — > IS transition. Using as a benchmark the value of the 
same amplitude for the 25 — > IP transition from the known rate of the decay ip' — > h c 7r° [1], 

1 The decays Xco — > r/ c ir + ir~ and Xd — > Vc^ are forbidden by parity. 



(1) 



r( Xc o^ c 7r°) = 3r(/i c ^ J/4>k°) ■ 



(2) 
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one can estimate T(xco — > Vc^ )) ~ (few keV), although the reliability of such estimate is 
not presently clear. 

The description within the QCD multipole expansion of the considered here pion transi- 
tions closely parallels that originally suggested [6] (also in the review [7]) for similar processes 
in bottomonium: Y(3>S) — > hb(lP)7T + 7T~ and T(3S) — > hb(lP)ir°, and an experimental ev- 
idence for the latter decay has been recently reported [8]. All these processes involve a 
transition with the change of both the total spin of the heavy hark-antiquark pair as well 
as its orbital momentum, so that they proceed through the interference of the El and Ml 
terms in the multipole expansion. In this respect these decays are quite different from the 
heavy quark spin-conserving processes, e.g. ip' — > J/i/jirir, or Xbji^P) Xfej(l-P) 7r7r ; an d 
using either of the latter for normalization [9] may lead to erroneous predictions. 

All essential details of the calculation can be found in the literature, e.g. in the review 
[7J. I briefly outline the reasoning here in order to make the discussion somewhat more self 
contained. The El and Ml terms in the multipole expansion are described by the following 
terms in the effective Hamiltonian 

H E1 = -lcr-E\ H M i = -7^— C (A • B a ) , (3) 
2 2m c 

where £ a = — ^2 * s the difference of the color generators acting on the quark and antiquark 
(e.g. t\ = X a /2 with A a being the Gell-Mann matrices), r is the vector for relative position of 
the quark and the antiquark, A = (cr c — (?c)/2 is the difference of the spin operators for the 
the quark and antiquark. Finally, E a and B a are the chromoelectric and chromomagnetic 
components of the gluon field strength tensor. The assumed here normalization convention 
is that the QCD coupling g is absorbed into the definition of the gluon field strength. 

The general expression for the amplitudes of the transition between the states X and Y 
of charmonium with emission of one or two pions can thus be written in a form factorized 
into the product of the quarkonium transition amplitude and the amplitude for production 
of the pions by the gluonic operator 

A(X Yn°) = -L- (7r°i^|0) (Y^inQ A k + A fe £r^ fe |X> , 

A(X Ftt+tt-) = ^- (7r+7r-|£?fl2|0) {Y\i\ Tl Q A k + A k g ri )?\X) , (4) 

where Q is the Green's function for propagation of the heavy quark pair in a color octet 
state. 
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The amplitude for the single tc° production by the gluonic operator in Eq.Q is deter- 
mined [TU] by the chiral anomaly in QCD once the isospin breaking by the mass difference 
Trid — m u between the down and up quarks is taken into account: 

/ OlnaDaln\ Z 2y/2 TT 2 m d - m u 2 

(n \ E i B k\ Q ) = 6 ik — ~ U ml , (5) 

3 m d + m u 

where f n ~ 130 MeV is the pion decay constant. The amplitude for the dipion production 
by the gluonic operator is evaluated by noticing that the the product E^Bf. is proportional 
to the off-diagonal (0m) component of the gluon field energy-momentum tensor 9^ v . The 
matrix element (v"7r|6^|0) has been studied by Novikov and Shifman [H] in terms of the 
fraction pa of the pion momentum carried by gluons, and the relevant amplitude in Eq.(j4]) 
is found as 

(7r+7r|£f5£|0) = 2na s p G e ikm (E lP2m + E 2 p lm ) , (6) 

where E\ , E 2 and p\ , p 2 are the energies and the momenta of the two pions. The estimated pTj 
numerical value of the coefficient in this expression, 2ira s pG ~ 2 4- 2.5 is in agreement with 
the measured behavior of sub-leading details (in particular the .D-wave) in the transition 
ip' — > J /ip-Kit (a discussion and further references can be found in Ref. [7]). 

The spatial structure of the matrix elements (jSJ) and (E]) uniquely determines the selection 
rules for the heavy quarkonium amplitudes in the expressions (J3J). Namely, the 8n~ dependence 
of the matrix element in Eq.© selects only the one-pion transitions between quarkonium 
states with the same J: AJ = 0, while the antisymmetric in % and k dependence eikm in 
Eq.(|6]) selects only the two-pion transitions with |AJ| < 1. For this reason neither the decay 
Xc2 Vcir , nor Xc2 Vc^ir receives any contribution from the leading El — Ml term in 
the multipole expansion, so that both these decays should be significantly suppressed. 

The amplitudes of all the allowed by the selection rules IP — > IS transitions can be 
written in terms of one radial overlap amplitude 

hs,ip = ttt~ — (Ris\rQp + Qsr\Rip) , (7) 
lo m c 

where R\s, R\p and Qs,p are the radial wave functions for quarkonium and the radial partial 
wave color-octet Green's function for the corresponding orbital momentum. The expressions 
for the transition amplitudes read as 

A[h c (a h ) -> J/i[)(a^) vr°] = 2y ^ 7T md mu f n m ^ j isip . ^) ; 

6 m d + m u 
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At v 0\ 2 y^7T 2 m d -m u 2 
-4(Xco ^VcK ) = 7 = ■ U m n I 1S1P , 

Alxdia^ ^7 lc 7i + 7i-} = 2V2XI ls , 1 p [E 1 (p 2 -a 1 ) + E 2 (p 1 -a 1 )] , (8) 

where a^, and c?i are the polarization amplitudes of the spin one particles hj,, J/ip and 
Xd respectively, and the notation A = na s pG is introduced. Numerically A ~ 1. 

Clearly, the unknown charmonium overlap integral I±s,ip cancels in the ratio of the rates 
of the decays listed in Eq.flHJ), and one readily finds the relation fl2]) between the first two 
decay rates (the phase space in these two decays is essentially the same), and 

r(Xc0^^ c 7r°) 4tt 2 ( 2 m d -m u 2 \ 2 13.7 

' 7T ■ f n m n p n o « — — , (9) 



r(Xci ->■ ^ c vr+7r-) A$ V m d + m u 11 V 17 A 2 

where p^o w 413 MeV is the pion momentum in the transition Xco — > Vc ^ and $ is the phase 
space integral for the decay Xd — > r] c 7i + n^: 

$ = / [EIpI + El pi) PlP2 dE 1 « (249 MeV) 7 , (10) 

with A = M(xd) — M(t] c ). The numerical value (m<j — m u ) / (m d + m u ) ~ 0.3 [12] is used in 
the estimate in Eq.flHJ). Naturally, at A ~ 1 this estimate gives the relation in Eq.flT]). 

It can be mentioned that the rate of one more process described by the same overlap 
integral hs,ip, the decay h c — > J/ipim, is strongly kinematically suppressed due to the 
significantly smaller mass difference between h c and J/ip. Indeed, within the same approach 
the amplitude of this decay is given by 

A[h c (a h ) -> J/i[)(af) tt+tt - ] = 2 A hs,xp e ikm {Ex p 2m + E 2 p lm ) a hi a^ h , (H) 

and the rate can be related to that of Xd ~~ Vc vr + 7r _ as 

where $i is the phase space integral as in Eq. lJTU]) . but with A = M(h c ) — M(J/ip). Numer- 
ically $i ~ (181 MeV) 7 . The value of the ratio ffl2|) and the relations (Q and (J2J) imply, in 
particular, that the rate of the transition h c — > J/ip7T + iT~ is approximately 40 times smaller 
than that for h c — > J/ipn in full agreement with the analysis of Ref. [6] of the transitions 
involving an 1 P\ quarkonium. 

The absolute scale for the rates of the discussed decays depends on the unknown overlap 
integral I\sxp- It might be helpful, for an approximate orientation, to compare these decays 
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with the transition if>' — > h c 7r°, whose rate T 0.26 keV [I] is determined by a similar 
overlap integral I\p t 2s- One thus readily finds the ratio of the rates in terms of the ratio of 
the corresponding /: 

r(Xco Vcn°) =3T(h c -> J/^tt ) = r(^ -»• M°) 7^ 

(13) 

where fc^o w 85MeV is the pion momentum in the decay ip' — > h c 7T°. Thus, if the overlap 
integrals for the 2S — > IP and IP — > IS transitions are of similar value, the rate r(Xco — > 
i] c tt°) should be in the ballpark of few keV, and the branching fractions in Eq.([T]) should be 
of the order of (few)xlO -4 . 

As uncertain as the present estimate of the branching fractions ([1]) is, it still indicates that 
the recent experimental upper bound pQ B(Xci Vc' n ~ + 'n~~) < 0.32% is by about an order 
of magnitude higher than a reasonable theoretical expectation. Also, the present estimates 
imply that it may be more advantageous for experimental studies to search for the decays 
with single pion: Xco ~^ Vc 71 ' and h c — > J/tpn . In spite of being suppressed by the small 
breaking of isospin by the light quark masses, these decays are enhanced by the contribution 
of the chiral anomaly (Eq.Q) as has been pointed out long ago [6] and as it is indicated 
by the E760 results [5] and by the data [8] on the transitions from T(3S) to hb(lP) in 
bottomonium. 

This work is supported, in part, by the DOE grant DE-FG02-94ER40823. 
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